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Abstract 

The present study of small molecules containing silicon has been motivated by (a) 
the considerable interest being shown currently in the kinetics and reactivity of such 
molecules, and (b) the biotechnological potential of silicon-derivate surfaces as substrates 
in the adsorption of, for instance, amino acids and proteins. Therefore, we have studied 
by (i) a semi-empirical approach and (ii) an ab initio procedure employing low-order 
M0ller-Plesset perturbation theory, the molecular correlation energies of some neutral 
closed and open shell silicon-containing molecules in the series SiX„Ym- Procedure (i) 
is shown to have particular merit for the correlation of the ionic members studied in the 
above series, while the ab initio procedures employed come into their own for neutral 
species. 
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1 Introduction 



In the last twenty years, much work has been carried out in order to calculate the 
correlation energy in atoms and molecules. In particular, in atoms, March and Wind 

showed how it is possible to recover the main trend of the semi-empirical correla- 
tion energies obtained long ago by Clementi J?, using simple arguments based on the 
correlation energy density (see also (3j). 

Within the framework of Density Functional Theory (DFT), many studies have been 
made to develop new correlation energy functionals. These are now available as a result 
of the work of Vosko, Wilk and Nusair (VWN) [3], Lee, Yang and Parr (LYP) [S], and 
Perdew and Wang (PW) 6 , and have all been introduced into some quantum chemistry 
packages (gaussian, nwchem, gamess, etc). 

In a previous study (hereafter denoted as Paper I) 0, some of the present authors 
have introduced a new semi-empirical technique in order to calculate the electron corre- 
lation energy in molecules. Essential to this work was the division of the total electron 
density of the molecule into parts belonging to specific nuclei, and a part belonging 
solely to the chemical bond and therefore, by using semi-empirical modelling, strictly 
correlated with the molecular bond-order (BO). One important conclusion of this work 
was that the molecular correlation energy is generally larger, in absolute value, than the 
sum of the correlation energies of the separated atoms. It is known from the pioneer- 
ing work of Pauling, Mulliken, Coulson and other authors [HI 1^ ITU], that bond order 
and molecular bond length are strongly correlated. In this sense our work is somewhat 
related to the work of Fulde et al. jlll I12j . where an analytical model of molecular 
correlation energies based on bond lengths was proposed. 

In Paper I we developed a parameterization for compounds containing only hydrogen, 
oxygen and carbon atoms and the related bonds X-Y (X, Y = C, H, O). The very good 
results obtained in the calculation of the correlation energy for some molecules containing 
these elements, compared with the experimental values, encouraged us to extend this 
approach to second-row elements and in particular to silicon-containing compounds, like 
fluorine and chlorine-substituted silylenes, silyl radicals and silanes. It is relevant in this 
context to note that Schlegel et al. I14| II 5j have calculated the heats of formation of 
these compounds by ah initio methods, introducing part of the correlation energy using 
the standard M0ller-Plesset perturbative expansion technique. 

We stress at this point that the kinetics and reactivity of small silicon-containing 
molecules are attracting considerable interest, due to the quite different reaction mech- 
anism of the silicon compounds with respect to the carbon compounds. For example, 
pyrolysis of silanes yields silylenes rather than silyl radicals, whereas photolysis produces 
silyl radicals jl6j . Moreover, silylenes and silyl radicals with various degrees of fluorina- 
tion (or chlorination) play an important role in chemical vapor deposition |17 | I18[ ITOj. 
But the main area of interest of these compounds are the silicon-derivate surfaces which 
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are widely used as substrates in the adsorption process of biological compounds, such 
as amino acids, proteins etc, in view of their application in different biotechnology areas 
such as biomaterials, biosensors, and bioseparation 20 . 

In this study we have investigated the correlation energies of SiX^Y^ (X, Y = H, 
F, CI; n + m = 4, 3, 2) compounds using both standard ab initio techniques, at various 
levels of approximation, and the bond-order method developed in Paper I. In section [21 
we describe briefly the theoretical methods used to estimate the molecular correlation 
energy. M0ller-Plesset theoretical techniques are summarized in subsection 12.11 while 
the Bond-Order Correlation Energy (BOCE) method is described in subsection 12.21 
Section |31 is concerned with the calculated correlation energies of SiX„Ym compounds 
and the ionization potentials (IP) of some of these molecules together with a comparison 
with the experimental values. Our final conclusions are recorded in Sectional 

2 Methods 

2.1 Ab initio methods 

Ab initio molecular orbital Hartree-Fock calculations (HF) were performed with the 
GAUSSIAN 03 package (2^ using the polarization 6-31G** basis sets. The 6-31G** basis 
contains a set of six d Cartesian functions on each of the heavy atoms, and a set of three 
p function on each hydrogen. Equilibrium geometries were obtained by full optimization 
using analytical gradient methods [221 ■ Electron correlation energy was estimated with 
M0ller-Plesset (MP) perturbation theory [2S1 up to fourth order (MP4), including all 
electrons and considering single, double excitations. Full molecular geometry optimiza- 
tion was performed up to third order M0ller-Plesset level (MP3). Quadruple excitations 
(MP4SQDT) were introduced by fixing the molecular geometry obtained from MPS 
calculations. For both HE and MP calculations, restricted Hartree-Fock (RHE) and un- 
restricted Hartree-Fock (UHF) methods were applied to closed and open shell systems, 
respectively. 

We stress, at the outset, that the emphasis of this work is the correlation energy. To 
present such results on a large number of Si-containing compounds, we have deliberately 
restricted the basis sets used. To extend these is straightforward for the future, should 
it be deemed instructive to do so. 

2.2 Bond-Order Correlation Energy (BOCE) Method 

The 'experimental' correlation energy Ef^P is defined as 

Ef^P = Es — EfiF, (1) 
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where Enp denotes the (molecular) Hartree-Fock energy and Es is the so-called Schrodin- 
ger (or "true") energy. In the case of atoms, the latter is obtained from the sum of the 
experimental ionization energies, plus Bethe-Salpeter and relativistic corrections. In the 
case of molecules, Es also includes smaller contributions, taking into account for the 
other molecular degrees of freedom, such as vibrational and rotational motions. In both 
cases, Es can be constructed from available data j24j following the standard procedure 
described in Cremer's papers [251 126j . By using this procedure and the available ex- 
perimental data, we have obtained Es for all molecules in the series SiX„Ym (X, Y = 
H, F, CI; n + m = 4, 3, 2). Although Ef^^ in Eq. is usually referred to as an 'ex- 
perimental' correlation energy, it obviously contains some theoretical input through the 
Hartree-Fock energy -EhFi whose degree of approximation can be controlled by choosing 
a sufficiently large basis set. Moreover, for some molecules, whenever the experimental 
data for the heat of formation, or the vibrational molecular frequencies were unavailable, 
we have made recourse to available theoretical estimates of these quantities, as indicated 
in the table captions. Fig. ^ summarizes the definitions of the experimental correlation 
energies, as well as of other various energies addressed in turn below. 

Within the bond-order approach [7j, the correlation energy can be estimated as 
follows. Given the formation reaction of a generic molecule, 

A + B — > C + binding energy, (2) 

the molecular energy can be partitioned as 

E{C) = E{A) + E{B) - binding energy = E{A) + E{B) + E{AB), (3) 

where E{AB) is the energy due to the A-B bond formation. Then, in general, the 
Schrodinger energy Es of a polyatomic molecule can be written as 

Es = J2Es{A)+ J2 EsiAB). (4) 

A all AB 

Likewise, in the HF approximation, the total molecular energy can be expressed in the 
form 

^hf = E^hf(^)+ E EMAB). (5) 

A all AB 

Subtracting Eq. Q from Eq. Q), one obtains the theoretical correlation energy as 
^theo ^ Y.iEs{A)-EnF{A)]+ [EsiAB) -EuFiAB)] 

A all AB 

= molecule) + E E'^'°{AB), (6) 

A all AB 

where we take as known the Schrodinger atomic energy from experimental data, as 
explained in the next section. In Eq. ©, -E^c^''(^in molecule) is the effective correlation 
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energy for the atom A in the molecule. This quantity is obtained from Ef^^^A) relative 
to an isolated atom A and rescaled taking into account the effective number of non- 
bonding electrons of the atom A in the molecule. In particular, if n is the total number 
of electrons of A involved in the bonds and Z its atomic number, -Ec^p (An molecule) is 
obtained as 

i^rP(^i„ molecule) = E^^^A)?^. (7) 

The next step in the present approximation scheme is to consider E^'^°{AB) as 
an analytical function of the bond-order Pab between the two atoms A and B of the 
molecule. In practice, we have assumed a linear form for E^'^°{AB) as a function of 
Pab, namely 

Ef'^^iAB) = aAsPAB, (8) 

where uab is a parameter which depends on the A~B bond. Moreover, the Pab bond- 
order, according to Lowdin's definition p7,, can be written as 

A B 

PAB = J2J2^f^-^t^-^ (9) 

with D = S^^'^TS^/'^ , S being the overlap matrix and T the first-order density matrix 
obtained from the HF calculation. One then obtains: 

^^'^^ = ^ i^r'' (An molecule) + E ^^AbPaB- (10) 
A all AB 

The parameters uab can be obtained from a model system which contains the A-B 
bond and for which is it possible to calculate the Schrodinger molecular energy Es and 
then the experimental correlation energy Ef^^. After obtaining the best values of the 
aAB parameters for each A-B bond, application of Eq. then permits to calculate 
the theoretical correlation energy E^^°. In Appendix^we report a specific example to 
illustrate the procedure used to obtain the uab parameters. 



3 Results and discussion 

3.1 Neutral closed and open shell SiX^Y„i compounds 

In order to calculate the correlation energy of fluorine and chlorine-substituted silylenes 
compounds, we have obtained the bond parameter uab of Eq. (jHI using the experi- 
mental correlation energy Ef^^ of some model molecules and the experimental atomic 
correlation energies. In Tab. ^ we report the Schrodinger and the HF atomic energies 
of the constituent atoms in their fundamental state and the related correlation energies. 



5 



In Tab. |21 we record results for some model systems together with the bond parameter 
extracted from Ef^P. All energy values are in Hartrees. 

In the second column of Tab. El we record the Schrodinger energy of the SiX^Y^ 
compounds under study, obtained from experimental values of the heat of formation 
and of the vibrational frequencies of these molecules jl31ll4l[T^ . For the compounds for 
which no such experimental data were available, the theoretical values of these quantities 
were used to calculate the Schrodinger energy. Columns 3-6 of Tab. 01 report the total 
energy obtained at the HF, MP2, MPS and MP4SDTQ4 (MP4) levels. In the last 
column, the BOCE energy is reported. For the sake of simplicity. Tab. O does not 
include the molecules employed in the BOCE treatment to calculate the bond-order 
parameter uab- Tab. 01 reports the experimental and the calculated correlation energies 
for the same molecules listed in Tab. |31 

From Tab.jni it is evident that using the ab initio methods, the successive corrections 
to the molecular energy are very small. The differences between the MP2 and MP3 
results range in general between 0.02 — 0.04 a.u., while the difference between the next 
successive corrections, MPS and MP4, is even smaller (0.01 — 0.02 a.u.). For some 
molecules (SiHFs, SiHF2, SiF2), the total energy increases on passing from the MP2 to 
the MPS correction. This trend is better evidenced by the correlation energies in Tab. 01 
There, the change in the values of Ec ranges between 0.02 and 0.05 a.u., on going from 
the MP2 to the MP4 correction. 

The very good agreement of the calculated molecular energy using the BOCE tech- 
nique with the Schrodinger value appears clearly from the values in Tabs. and ^ The 
highly accurate values of the correlation energy obtained from BOCE for all molecules 
of the series are confirmed from the calculated percentage error with respect to the ex- 
perimental correlation energy, ranging within 0.002 — 4.06 %. On the contrary, ab initio 
correlation energies yield very high percentage errors. In Tab.jSl for some molecules, the 
calculated correlation energy with the BOCE method gives a value higher than the ex- 
perimental one, and the corresponding molecular energy is lower than the Schrodinger 
energy. The absolute percentage error for these molecules has a value in the range 
between 0.002 % and 0.4 %. However, it is important to add some comments here. 
The first one is that the Schrodinger energy is obtained from experimental data and 
any experimental measurement is subject to an error that, in general, is about ±5 %. 
The second comment, pertaining to the compounds corresponding to a molecular energy 
lower than the Schrodinger energy, is that one needs to note that in closed shell systems, 
such as those considered here, this happens only in two cases, where we have used some 
theoretical values (vibrational frequencies, or molecular heat of formation, or both) in 
order to determine the Schrodinger molecular energy. 

The results of Tab. jH are shown in Fig. [21 From this figure it can be seen that 
the trends of the calculated correlation energies using ab initio methods are close to 
the experimental behaviour, but the values are on average 1.5 a.u. higher than the 
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experimental values. Moreover, the increase in the calculated correlation energy on 
passing from MP2 to MP4 is negligible and, taking into account that the computational 
effort increases considerably from MP2 to MP4, we can conclude that it is not necessary 
to make corrections at higher order than MP2. This is confirmed by Fig. |31 where the 
absolute percentage error between the experimental and theoretical correlation energies 
is almost constant for each compound on passing from MP2 to MP4. 

The extent of the agreement between the experimental and the calculated BOCE 
correlation energies is shown in Fig. [21 (last panel), where the experimental and the 
BOCE curves are almost superimposed and hardly distinguishable. Consequently, as 
shown in Fig. 01 the percentage error between the experimental and the calculated 
BOCE correlation energies is nearly constant, and varies between 0.002 % and 4.0 %. 

3.2 Ionization potentials (IP) 

We have also calculated the ionization potentials (IP) of some SiX„Ym compounds, 
using both ab initio and BOCE methods. These results have then been compared with 
the experimental values. Within the so-called ASCF procedure, the theoretical values 
of IP (in eV) are defined as the difference between the total energies of the neutral and 
the ionized molecule, i.e. 

IP = E-E^ (11) 

(see also Fig. Table 13 reports the experimental and calculated IP for the molecules 
under study. In these calculations, the geometries of both the neutral and the ionic 
species have been optimized. 

In the second column of Tab. Owe record the IP values obtained as the difference E — 
E~^ , both calculated at HE level (without correlation), while in columns 5, 7, and 9 this 
difference refers to the calculated values at the MP2, MPS, and MP4 levels, respectively. 
Finally, column 11 reports the calculated IP values using the BOCE method. Fig. |3 
displays the data in Tab. O 

From Tab. El it is clear that the MP corrections, contrary to the case of the neutral 
compounds, are important in the IP calculation. In fact, in all compounds of this series 
(with the exception of SiH4), the absolute percentage error of the IP at the HE level is 
very high (~ 7 %; see Fig. jSJ, and the introduction of correlation effects at MP2 level 
reduces this error to 3.26 %, on the average. Then the introduction of correlation in the 
estimate of IP yields a calculated value closer to the experimental one. The trend of 
the next corrections, both at the MPS and at the MP4 levels, is not constant. For some 
compounds of this series, SiH4, SiIl2F2, SiHCls, and SiF4, the percentage error decreases 
on passing from MP2 to MP4. On the contrary, for other compounds (SiCl4, Sills CI) 
the error increases, while for SiHFs and SiIl2Cl2 the best estimate of IP is obtained at 
the MPS level. The average absolute percentage errors for the MP2, MPS, and MP4 
corrections are S.26 %, S.26 %, and S.44 %, respectively, thus showing that all these 
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low order MP corrections give almost equivalent approximations of the IP. Coming now 
to the BOCE approach, from the absolute percentage errors plotted in Fig. [S] one may 
conclude that this method yields more accurate estimates of the IP, on the average, than 
the ab initio methods, although the MP methods occasionally produce lower absolute 
percentage errors in IP than the BOCE method. 

A final consideration concerns the calculated differences. A, between the experimen- 
tal and calculated IP, reported in Tab. [51 Fig. ^schematically defines the experimental 
and theoretical correlation energies for both neutral and ionic molecules. According to 
Fig. ^ we may then deduce that 

jpexp _ jptheo ^ (^exp _ ^theo^ _ (^+cxp _ ^+thco^ = A - A+, (12) 

where a superscript + refers to the same quantity in the ionic molecule. In Tab. El 
we record the experimental and theoretical ionization potentials (IP), along with their 
difference, A. The fact that the average error in the calculation of the IP is not very 
different within the MP and BOCE methods, is an indication that there are large can- 
cellations between the correlation energies of the neutral and the ionic molecules. In 
turn, this means that within the BOCE approach also the correlation energy of the ions 
is well approximated by using the same bond parameters a^s, as derived for the same 
model molecules. 



4 Summary and concluding remarks 

The biotechnological possibilities of silicon-containing molecules in substrate layers, plus 
much current interest in understanding the quite different reaction mechanism of silicon 
compounds in comparison with carbon compounds, has motivated us to extend our 
earlier work [71 on the latter class to the case of Si. 

What seemed to us important in this different series of small molecules was to assess 
the utility of our earlier semi-empirical use of molecular bond-order to estimate elec- 
tron correlation energies in some members of C-containing molecules when applied to 
SiX„Ym, where X, Y = H, F, CI, and n + m = 4, 3, 2. In this process, we have found it 
valuable to compare and contrast our bond-order approach with some ab initio results 
we have obtained using low-order M0ller-Plesset perturbation theory. While, by these 
two approaches, we record quite a number of useful results for molecular correlation 
energies of Si-containing molecules, we wish especially to stress that the semi-empirical 
bond-order approach proposed in [3| and applied there to C-containing molecules con- 
tinues to be valuable for molecules of the class SiX„Ym. It seems particularly useful 
when ionicity plays a role, whereas the low-order M0ller-Plesset approach comes into its 
own for essentially neutral members of this series. 

In conclusion, we emphasize the essence of the present approach which is to yield an 
impressive empirical correlation for obtaining bond additivity corrections to the energy 
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based on bond orders calculated from a Lowdin population analysis. In order to present 
rather extensive results for a series of Si-containing molecules, we have, in the present 
work, accepted the limitations of the small basis set used. However, should it prove 
instructive for future purposes, it is a straightforward matter to extend the basis, even 
if somewhat time-consuming. Plainly, then, some quantitative changes will occur and 
also some specific technical points can be examined such as the numerical modifications 
in the Lowdin population analyses, and the parameters obtained from them, due to 
changes in the basis set. Another matter deserving attention if eventually larger basis 
sets are used concerns specifically the Si-F bond, which present indications suggest is 
somewhat problematic in the calculations, or experiment, or perhaps both. 

We are encouraged by the results of the present investigation to attempt a pilot 
study of the relevance of our bond-order considerations to the still more difficult area of 
the kinetics and reactivity of small Si-containing molecules. We hope to report on this 
area, approached via our semi-empirical procedures, at a later date. 



A Calculation of the bond-order parameter asin 

In this Appendix we describe the procedure employed to estimate the bond parameter 
aAB in Eq. 0, when AB = SiH. According to Eqs. ^ and (jH]), the experimental 
correlation energy of the SiH4 molecule can be written as 

ET^iSm,) = E^iSkn molecule) + E ET^iB^^^^^J + agiH E P^, + a^U E P^f^, 

j=l 1=1 1=1 

(13) 

where PsiH and Phh and the bond-orders for each Si-H and H-H bonds in SiH4, respec- 
tively. From the II2 molecule one immediately obtains the bond parameter for the H-H 
bond as 

E^^P(R2) — 2ii'c''^(Hin molecule) ,x 

aun = . (14) 

From Eq. with Z = n = 1, one has -E'^^P(Hin molecule) = 0. Since in the hydrogen 
molecule Phh = 1, one finds 

aHH = £^c'P(H2) = 4.28 • 10-2. (15) 

In Tab. iniwe report the symmetric matrix of the bond-order for SiH4. The diagonal 
elements are the sum over the off-diagonal elements, and represent the total electrons 
involved for each atom in the bonds. 

Using the experimental value -E^^P(SiH4) = 2.1742 a.u., and employing the results 
of Tables ^ and IHl from Eq. (fT^ we eventually find 

asiH = 6.238 • 10-^ (16) 
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as quoted by Tab.HJ 
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Atom 


S 




Es 


rncxp 




Si 


3 


-2.888318E+02 


-2.893116E+02 


4.798144E- 


-01 


H 


2 


-4.982329E-01 


-4.993000E-01 


1.067100E- 


-03 


F 


2 


-9.936496E+01 


-9.972500E+01 


3.600440E- 


-01 


CI 


2 


-4.594480E+02 


-4.600301E+02 


5.821385E- 


-01 



Table 1: Hartree-Fock energy, £^hf, estimated with the 6-3 IG** basis set, Schrodinger energy, 
Esi and experimental correlation energy, i^^^P, for Si, H, F and CI atoms. The second column 
lists the values of the spin multiplicity, S. All values are in a.u. 



A- 


-B 






H 


H 


H2 


4.287390E- 


-02 


Si- 


H 


SiH4 


6.237950E- 


-02 


F- 


-F 


F2 


1.674987E- 


-01 


Cl- 


-CI 


CI2 


1.314397E- 


-01 


H 


-F 


HF 


1.137677E- 


-01 


H 


CI 


HCl 


7.700011E- 


-02 


Si- 


-F 


SiF4 


1.273342E- 


-01 


Si- 


CI 


SiCU 


1.024566E- 


-01 


F- 


CI 


FCl 


1.586152E- 


-01 



Table 2: Lowdin bond parameters aAB for A-B bonds in several ABn model molecules. See 
Appendix El for the derivation of asiH- For the FCl molecule (last row), the experimental 
vibrational frequency being unavailable, a theoretical value has been used in the calculation 
of ^5. 
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Molecule 


Es 


E-RF 


MP2 


MP3 


MP4 


BOCE 


SiHgF 


-391.149094 


-390.152840 


-390.438375 


-390.450934 


-390.463879 


-391.136150 




-490.474680 


-489.084835 


-489.538902 


-489.543176 


-489.563139 


-490.460715 


SiHFg 


-589.799528 


-588.019851 


-588.642790 


-588.639082 


-588.665926 


-589.792771 


SiHgCl 


-751.377125 


-750.187745 


-750.434135 


-750.461123 


-750.469475 


-751.374545 


SiHsCls 


-1210.934066 


-1209.143833 


-1209.520705 


-1209.553336 


-1209.564921 


-1210.929943 


SiHCla 


-1670.491322 


-1668.097366 


-1668.607261 


-1668.644964 


-1668.660441 


-1670.486756 


SiHaFCl"'^ 


-850.697705 


-849.113982 


-849.529345 


-849.547673 


-849.563687 


-850.694892 


SiHFsCl'^'* 


-950.024261 


-948.045333 


-948.630412 


-948.640346 


-948.663824 


-950.024303 


SiHFCb"'^ 


-1310.255828 


-1308.071123 


-1308.610000 


-1308.642327 


-1308.662046 


-1310.255383 


SiFsCls'^''' 


-1409.578850 


-1406.998802 


-1407.716386 


-1407.731725 


-1407.759092 


-1409.583792 


SiFgCl" 


-1049.355060 


-1046.974456 


-1047.728506 


-1047.730616 


-1047.761047 


-1049.353836 


SiFClf 


-1769.814012 


-1767.022899 


-1767.704115 


-1767.732799 


-1767.756713 


-1769.813079 


SiH2 


-290.561953 


-290.002631 


-290.093976 


-290.111950 


-290.117816 


-290.539281 


SiHF 


-389.887310 


-388.933412 


-389.140967 


-389.205213 


-389.219847 


-389.884142 


SiFa 


-489.218209 


-487.884672 


-488.321054 


-488.320264 


-488.343220 


-489.206212 


SiHCl^ 


-750.122200 


-748.9442 


-749.195368 


-749.219098 


-749.228569 


-750.077830 


SiCl2 


-1209.699563 


-1207.943683 


-1208.301785 


-1208.330378 


-1208.344114 


-1209.678803 


SiH| 


-291.173027 


-290.610579 


-290.709229 


-290.726369 


-290.731504 


-291.175489 


SiH2F« 


-390.493513 


-389.525900 


-389.791422 


-389.800054 


-389.812671 


-390.480470 


SiHFf 


-489.813614 


-488.452373 


-488.886854 


-488.886530 


-488.906703 


-489.800833 


SiHsCF'* 


-750.723808 


-749.566110 


-749.792607 


-749.815941 


-749.823949 


-750.688190 


SiHClg'^ 


-1210.279079 


-1208.521554 


-1208.878868 


-1208.907693 


-1208.919212 


-1210.281674 


SiFa 


-589.132214 


-587.381225 


-587.984340 


-587.975648 


-588.002862 


-589.126715 


SiClf 


-1669.840418 


-1667.474563 


-1667.964920 


-1667.998734 


-1668.014351 


-1669.838233 



Table 3: Schrodinger and theoretical energies (in a.u.) for open and closed shell of SiX„Ym 
compounds. The upper table refers to relevant energies of closed-shell silicon compounds: 
especially the bond-order correlation energy. The lower table refers to the energies of open- 
shell silicon compounds. Notes : (a) The Schrodinger energies have been calculated using 
theoretical values of the vibrational frequencies, (b) The Schrodinger energies have been 
calculated using theoretical values of the molecular heat of formation. 
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ivioiecuie 


pexp 




i\/rp^ 


lVlr^4 




oiJi3r 


u.yyozoo 


U.ZoOOoO 


n OQsno/i 
u./youy4 


u.oiuyoy 


u.yoooiu 


&lJl2-t 2 


i.ooyo4D 


U.404UDO 


U.4ooo4i 


U.4 ( oooO 


i.oi OooU 




i. / l\)0IO 


U.ozzyoy 


u.oiyzoi 


U.040U 1 


1. 1 1 zyzU 


SiH'^Cl 


1 189380 


246389 


273378 


281730 


1 186800 




1.790233 


0.376872 


0.409502 


0.421088 


1.786110 


SiHCla 


2.393956 


0.509894 


0.547597 


0.563074 


2.389390 


SiH2FCF''' 


1.583723 


0.415363 


0.433691 


0.449706 


1.580910 


SiHF2CP''' 


1.978928 


0.585079 


0.595013 


0.618492 


1.978970 


SiHFCl2''' 


2.184705 


0.538877 


0.571204 


0.590923 


2.184260 


SiFof!!"''' 


9 580048 


71 7584 


739Q94 


760991 


9 584990 


SiFaCP 


2.380604 


0.754049 


0.756160 


0.786592 


2.379380 


SiFClg 


2.791113 


0.681216 


0.709900 


0.733814 


2.790180 


SiH2 


0.559322 


0.091345 


0.109319 


0.115185 


0.536650 


SiHF 


0.953899 


0.207555 


0.271801 


0.286436 


0.950730 


SiF2 


1.333537 


0.436382 


0.435591 


0.458548 


1.321540 


SiHCl'' 


1.177958 


0.251168 


0.274897 


0.284368 


1.133630 


SiCl2 


1.755880 


0.358102 


0.386695 


0.400432 


1.735120 


SiHf 


0.562448 


0.098650 


0.115791 


0.120925 


0.564910 


SiH2F'^ 


0.967613 


0.265521 


0.274154 


0.286771 


0.954570 


SiHFf 


1.361241 


0.434481 


0.434157 


0.454329 


1.348460 


SiH2Cl'^'* 


1.157698 


0.226497 


0.249830 


0.257838 


1.122080 


SiHClg''' 


1.757525 


0.357315 


0.386139 


0.397659 


1.760120 


SiFs 


1.750989 


0.603115 


0.594423 


0.621637 


1.745490 


SiClf 


2.365855 


0.490357 


0.524171 


0.539788 


2.363670 



Table 4: Experimental and theoretical correlation energies (in a.u.) for open and closed shell 

of SiX„Y„j compounds. The upper table refers to closed-shell, whereas the lower table refers 
to open-shell compounds. Notes : (a) The Schrodinger energies have been calculated using 
theoretical values of the vibrational frequencies, (b) The Schrodinger energies have been 
calculated using theoretical values of the molecular heat of formation. 
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Molecule 


jpexp 


HF 


A 


MP2 


A 


MP3 


A 


MP4 


A 


BOCE 


A 


SiH4 


11.60 


12.38 


-0.78 


10.70 


0.90 


10.75 


0.85 


10.78 


0.82 


11.16 


0.44 


SiHsF 


11.70 


12.20 


-0.5 


11.52 


0.18 


11.67 


0.03 


11.99 


-0.29 


11.49 


0.21 


SiHgCl 


11.40 


10.37 


1.03 


11.40 


0.00 


11.08 


0.32 


11.08 


0.32 


11.53 


-0.13 


SiHaFs 


12.20 


11.77 


0.43 


12.07 


0.13 


12.24 


-0.04 


12.11 


0.09 


12.45 


-0.25 


SiHsCb 


11.70 


11.36 


0.34 


11.52 


0.18 


11.56 


0.14 


11.43 


0.27 


11.86 


-0.16 


SiHFs 


14.00 


12.43 


1.57 


12.78 


1.22 


13.02 


0.98 


12.88 


1.12 


13.51 


0.49 


SiHCla 


11.40 


12.04 


-0.64 


11.67 


-0.27 


11.68 


-0.28 


11.56 


-0.16 


12.02 


-0.62 


SiF4 


15.70 


17.82 


-2.12 


15.12 


0.58 


14.93 


0.77 


15.35 


0.35 


16.54 


-0.84 


SiCl4 


11.80 


12.53 


-0.73 


11.48 


0.32 


11.44 


0.36 


11.34 


0.46 


11.85 


-0.05 



Table 5: Experimental and theoretical ionization potentials (in eV). The A columns report 
the differences between the experimental and the theoretical values. 





Si 


H 


H 


H 


H 


Si 


3.970845 










H 


0.993248 


1.015563 








H 


0.993248 


0.006901 


1.015563 






H 


0.991655 


0.007873 


0.007873 


1.014972 




H 


0.992695 


0.007541 


0.007541 


0.007571 


1.015348 



Table 6: Bond-order matrix of SiH4. 
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Figure 1: Schematic determination of the experimental and theoretical correlation energies 
for both neutral and ionic molecules. 



16 




5 10 15 20 5 10 15 20 

Figure 2: Correlation energies Ec (in a.u.) of closed and open shell silicon compounds, for 
the 24 molecules in Tab. EJ The abscissa is the row index in Tab. |3 The vertical dashed line 
separates closed from open shell molecules. In all insets, open squares refer to experimental 
values, while open circles refer to calculated values. 
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Figure 3: Absolute percentage errors, 100|(E;:^p - E*^"°)/E;:^p|, of the theoretical correlation 
energies with respect to the experimental correlation energy in Tab.|3]and Fig.|21 The abscissa 
is the row index in Tab. 01 The vertical dashed line separates closed from open shell molecules. 
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Figure 4: Ionization potentials IP (in eV) for the 9 Si-containing molecules in Tab. El The 
abscissa is the row index in Tab. |31 Open squares refer to experimental values, while open 
circles refer to calculated values. 
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Figure 5: Absolute percentage errors, 100|(IP""P - ipth<=o)/ip<=-P|, of the theoretical ionization 
potentials with respect to the experimental values in Tab. and Fig. EJ The abscissa is the 
row index in Tab. El 
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